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The self-gravitating gas in thermal equilibrium is studied using a Newtonian potential regularized 
at short distances. This short distance cutoff permits us to obtain a complete description of the 
gas including its collapsed phase. We give a field theory description of the iV-body regularized self- 
gravitating gas in the canonical ensemble. The corresponding functional integral is dominated in the 
N — > oo limit by saddle points which provide a mean field description. The well-known dilute solu- 
tions (isothermal spheres) are recovered. We find new solutions which are regular in the regularized 
theory but become singular in the zero cutoff limit. They describe collapsed configurations where 
the particles are densely concentrated in a region of the size of the cutoff. These collapsed solutions 
provide the absolute minimum of the free energy. We find further new solutions which interpolate 
between the collapsed and the dilute configurations and describe tunneling processes where the gas 
collapses. The transition probability for such collapse processes turns out to be extremely small for 
large N. That is, the dilute solutions are in practice stable in the regime where they are locally 
stable. 



Contents 

I. Introduction Q] 

II. Mean Field Approach to the Self- Gravitating Gas 

III. Spherically symmetric solutions 

IV. Solutions Regular as A — > @ 



10 
13 
15 



V. Solutions Singular as A — > 0. 

A. Collapsed phase 

B. Interpolating solution: metastability 
References [l? 

I. INTRODUCTION 

Selfgravitating system are ubiquous in the universe. Selfgravity governs the formation and evolution of the large 
scale structure in the universe as well as the dynamics of galaxy formation, stars and interstellar medium [3, ■ It 
is therefore of great interest to study self-gravitating fluids in thermal equilibrium although most astrophysical and 
cosmological systems are not in exact thermal equilibrium [H-Q- 

The self-gravitating gas has peculiar properties from the point of view of equilibrium statistical mechanics: the 
most important configurations are inhomogeneous and the thermodynamic functions exist in the dilute limit [H, B 0] 

N 

N^^, V^w, fixed, (1.1) 
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where N is the number of particles and V stands for the volume of the box containing them. In such a limit, the 
internal energy, the free energy and the entropy turns to be extensive. That is, they take the form of N times a 
function of intensive dimensionless variables. In the case of a spherical box of radius R (so that V = 4irR 3 /3) there 
is only one such variable, that is the ratio of the characteristic gravitational energy Gm 2 N/R (m is the particle mass 
and G is Newton's constant) and the kinetic energy of the order of the temperature T of a particle in the gas [EIE01 : 

Gm 2 N 

For more complex geometries there will be, besides rj (in such cases R is to be identified with the linear scale of the 
box), other dimensionless shape parameters. 

In this paper we restrict the analysis to the simplest case of a spherical box and follow the approach of refs.^,0,0] 
where the starting point is the partition function for non-relativistic particles interacting through their gravitational 
attraction in thermal equilibrium. The configurational part of this partition function is expressed as a functional 
integral over a scalar field 4>{x) [proportional to the the gravitational potential] with an Euclidean action (that is 
minus the logarithm of statistical weight) proportional to N . Therefore, the saddle point technique (or mean field 
method) permits to compute the configuration partition function in the large N limit. 

As investigated in refs.Jl, d, @, 0] the self-gravitating gas in the canonical ensemble is in a gaseousphase for 
< rj < t\t where rjr = 2.43450. . .. (Further studies on self-gravitating gases are reported in ref. 0, [1, [t| (HI). 
At this point the isothermal compressibility diverges and the speed of sound inside the sphere becomes imaginary. 
This clearly shows that the gas phase collapses at this point in the canonical ensemble as a consequence of the Jeans 
instability. Monte Carlo simulations explicitly show that the gas collapses for rj > r\T into a very dense and compact 
phase Q. 

It must be recalled that the small fluctuations around the mean field solution in the canonical ensemble are stable 
for 77 < r/c = 2.517551 . . .@. Namely, fluctuations in the partition function do not feel the divergence of the isothermal 
compressibility neither the fact that the speed of sound becomes imaginary. 

Even in the regime where the gaseous phase is locally stable there are always collapsed configurations which have 
much lower energy than the gas. In fact, the energy has its absolute minimum in the collapsed phase where all particles 
fall to one point and this configuration completely dominates the free energy, causing its divergence to minus infinity. 
Of course, interactions other than gravitational always dominate for short distances. Therefore, the properties the 
collapsed phase depend on the short distance physics. This physics is different according to the nature of the 'particles' 
considered. They can be from galaxies till molecules or atoms. Still, we can study the regime when collapse happens 
just introducing a short distance cutoff (ultraviolet regulator) to the gravitational interaction. 

We present a complete description of the self-gravitating gas, including the collapsed phase, using a Newtonian 
potential regularized at short distances. We choose the Newton potential to be — G m 2 /A for interparticle distances 
r < A instead of Newton's law — G m 2 /r. The cutoff A is chosen very small compared with the size of the box. 

We give an exact mapping of the iV-body regularized problem into a field theory functional integral. This functional 
integral is dominated by its saddle points in the N — * 00, namely by the mean field equations. The mean field equations 
coincide for zero cutoff with the Lane-Emden equation obtained from hydrostatics and assuming an ideal gas equation 
of state. 

In this way we reproduce the well-known isothermal sphere solutions describing a self-gravitating gas in local thermal 
equilibrium. These diluted solutions get small corrections for small cutoff. 

We find a host of new solutions (saddle points) from our regularized mean field equations. These new solutions are 
singular in the limit of zero cutoff and cannot be find from the standard (no cutoff) Lane-Emden equation. In these 
new solutions, which are perfectly regular for non-zero cutoff, a large number of particles concentrate near the origin. 
There are two new types of solutions. 

The first type are the collapsed saddle point solutions: all the particles are densely concentrated in a region of size 
of the order of the cutoff. These solutions are studied analytically and provide the absolute minima of the free energy, 
which reads, at dominant order for small cutoff A <C R: 

Gm 2 N(N-l) , R 3 

F - F a ~ \ '- + NT log ■ 



2 A a {A/2) 3 ' 

Here, T is the temperature and Fq is the free energy of an ideal gas. The first term is just the potential energy of N 
particles clustered in a small sphere of radius A/2, where the regularized gravitational interaction is the same for all 
^N(N — 1) particle pairs. The second term is T times the entropy loss in the collapse. This free energy is large and 
negative, unbounded from below as A — > 0. In particular, this free energy is well below the free energy of the dilute 
solution. 

The second type of saddle point solutions interpolate between the gas phase and the collapsed phase, in a sort of 
pre-collapse. This solution has a finite action in the zero cutoff limit and describes the tunneling transition where a 
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small fraction of particles coalesce into a region of the size of the short-distance regulator so that the density near 
the origin is very large. The tunneling probability to collapse for a dilute solution turns out to be extremely small for 
large N and small A. This means in practice that dilute solutions which are locally stable can be considered stable. 
As mentioned above this happens in the canonical ensemble for < r\ < t\t where t)t — 2.43450 . . .. 

The parameter rj for the collapsed solutions can be arbitrarily large. That is, the particle density does not need to 
be dilute N ~ [see eq. (|l.l|) ] but we can have N ~ V. 

We obtain in sec. IVBI an estimate for the lifetime for the dilute phase of the self-gravitating gas: 

r^-^log-e- ^ ~At 0g A e ^l2G^r g A)\- 

One can see that the lifetime becomes infinitely long in the zero cutoff limit as well as when N — * oo at fixed cutoff 
[recall that R ~ N in the dilute limit of eq. 

In summary, we provide through mean field theory a complete statistical description of the self-gravitating gas 
including the absolute minimum of the free energy. 



II. MEAN FIELD APPROACH TO THE SELF-GRAVITATING GAS 



At short distances, the particle interaction for the self-gravitating gas in physical situations is not gravitational. 
Its exact nature depends on the problem under consideration (opacity limit, Van der Waals forces for molecules etc.). 
We shall just assume a null short -distance pair force; that is, we consider the Hamiltonian 



N 2 
PI 



^-E^-E 



Gm 2 



3 = 1 



j<k 



[ W ~ Qk\A 



where 



\Q\A = 



\q\ , for|g|>A 
A , for |q| < A 



(2.1) 



and A <C R is the short-distance cut-off. 

It must be stressed that the results presented in this work using the soft-core cut-off eq. (|2.1l) will be qualitatively 
the same for other types of cut-off since R 3> A. 

The partition function in the canonical ensemble reads 



Zn 



1 



d 3N P 



N\ h 3N J (2tt) 3N J v n 



where T = (3 1 is the temperature. Computing the integrals over the momenta pj, (1 < j ; < iV), 



d 3 p 

(2tt) 3 



exp 



2 m 



3/2 



yields 



N 



Zr f = 4r/ d™q TT expf, Gm ^ ) 



(2.2) 



(2.3) 



where 



7 (0) 

J N 



2nm 



3N/2 yN 



exp 



is the partition function of the ideal gas. Thus, 



F = F W _ T l g^ onf 



is the full canonical free energy. 



4 



We can now replace the integration over the particle coordinates with a functional integration over the configurations 
of the gravitational field ip produced by the particles. We start from the fundamental property of Gaussian functional 
integrals, namely 



Dip exp < — S[ip] + / d 3 x J(x) ip{x) > = exp 



~fd 3 xf d 3 x'J{x) g(x,x') J(x') 



(2.4) 



where S[ip] is the quadratic action functional 

1 



S[tp] = - / d 3 x / dV (p(x)g- 1 (x,x')(p{x') 



and the integration measure is assumed to be normalized such that J Dip exp{— <S[y>]} = 1. Then we set 

ni i\ G ft 1 
G{X,X = — 

\x - x'\ A 

and we identify the source field J(x) as (3m times the microscopic particle density 

N 

J(x) = /3mJ2 si3) ( x -Q]) ■ 

3=1 



(2.5) 



(2.6) 



Inserting eqs. (|2.5p and (|2.6p into eq. (12. 4p yields, 



N 



Dtp exp |^ — S[(p] + (3 m <p(Qj) J = exp 



L Gm 2 (3 
''^\lj-Qk\A 



Om- fi N f-f 

= e 23 II exp 

l<j<k<N 



Gm 2 (3 

\<lj ~ Qk\A 



(2.7) 



This provides a functional integral representation of the product over particle pairs in the integrand of eq. (|2.3p . 
Inserting the identity eq. (|2.7p in eq. (|2.3p and exchanging the functional integration with the integration over the 
coordinates qj we obtain 



*7Conf 
N 



e ™a J Dtp exp{-5[^]} J^^Q ex P 



N 



(2.8) 



where 



S[<p] 



{3 



8ttG 



d x <p V A ip 



and = -(4?rG//3) Q- 1 satisfies 



Since the qj are dummy variables, we have 



d 3N q exp 



\x - x'U 



N 



^5<- 3 H 



x — x 



Pm^2ip(qi) 



i=i 



and the configuration partition function becomes 



Z% ni = e I Vtpe 



d 3 f3m v (q) 



Sesl<p] 



N 



where 



S c g[p] = S[tp] + 5mt[v] 



SintM = —TV lOE 



d 3 x e ~ Pmv{x) 



(2.9) 



(2-10) 
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We want to stress that this field-theoretic functional formulation of the self-gravitating gas is fully equivalent, for any 
N, to the original one in terms of particles. The functional integral representation of the partition function presented 
in ref. for the canonical ensemble is slightly different but equivalent to eqs. (|2.9p and (|2.10p for large N. Let us 
also observe that the nonlocal operator V 2 A reduces to the standard Laplacian operator in the limit A — » 0. However, 
this limit is not so straightforward in the functional integral above, since the interaction term S- lri t in the action may 
(and indeed does) become unbounded from below in such limit. 
We now pass to dimensionless variables by setting r = x/R and 

—(3m <p(x) = 4>{r) . 

The action for 4>(x) now reads 



SoffM = Ns[4>] , a[4>] = -— I d'r <jy - log 



d 3 r 

r\<l 



log^, (2.11) 



where 



Gm 2 [3N A 
' R R 



-V2 1 - i -7T=^(r-r'). (2.12) 



We recall that the variable rj is the ratio of the characteristic gravitational energy Gm R N and the kinetic energy 
~ T = /3 _1 of a particle in the gas. For rj = the ideal gas is recovered. 

Having factored out of the action the number of particles N, it is natural to use the saddle points method to 
evaluate the functional integral over <fi in the limit N — > oo at fixed ?y, that is in the dilute limit. Strictly speaking, 
however, the action per particle s[4>] still depends on N through a = A/R oc A/N. Thus we should better regard the 
saddle point method as a mean field type approximation yielding the free energy F to leading order in N, as we shall 
now show. The stationarity condition 

_ n 

8<t>{r) 

leads to the regularized self-consistent Boltzmann-Poisson equation 

V^+^e* = , Q= [ d 3 re^ r) , for |r| < 1 , and 

Q J\r\<l 

V 2 J = , for |r|>l . (2.13) 

In the limit a — > this becomes the standard self-consistent Boltzmann-Poisson (or Lane-Emden) equation derivable 
from hydrostatic plus the assumption of a local, ideal gas equation of state [|| . 
We may rewrite eq. (|2.13l) in integral form as 

#r) = V [ d 3 r> , (2.14) 

J \r — r'\ a 

where 



p(r) = —— for \r\ < 1 , p(r) = for |r| > 1 , / p(r) = 1 , (2.15) 

Lis 

is the normalized particle density associated to (j)(r). 

We may also use the identity involving the standard Laplacian V 2 , 

-V 2 -^— = ±5(\r-r'\-a), 



to cast the field equation in the integro-differential form 

V 

Q a 



V 2 ^(r) = / d 3 R'S{\r-r'\ -a) e^ r ' ] . (2.16) 



k'l<i 
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Notice that the laplacian in the origin is always regular for nonzero cutoff 

H- f M\r'\ - n\ X 1 -') 



V 2 0(O) = --^/ dV 5{\r'\ - a) e^ r > , (2.17) 

Q a J\r'\<l 

implying, together with the fmitencss of 4>{0), that V(/> = in the origin. 

We can recast the action of the saddle point eq. (|2.11|) using the regularized equation of motion (|2.13|) as follows, 

»y>] = Tn I d 3 r0(r)e^)-log^ . (2.18) 

2 Q J\ r \<i 47T 

The free energy for large N can be written in terms of the stationary action as[f| 

F = F +N^-+NTs(r),a) + O(l) (at fixed a) , (2.19) 

where s(i],a) = s[4> s ], (f> s is a solution of eq. (|2.13|) , and F is the free energy for the ideal gas. Notice that, since a 
vanishes as ./V -1 as N — > oo at fixed A and r/, we can regard N T s(r],a) as extensive in the particle number, that is 
linear in N, only if it is regular as a — + 0. This holds true for stationary point which are regular at the a = limit. We 
show below that there are saddle points which become singular in the a — > limit and consequently the free energy 
is not proportional to N . Similarly, the subleading terms in eq. (|2.19|) are order one only for saddle points which are 
regular for a = 0. What is always true is that they are indeed subleading for any given saddle point. 



III. SPHERICALLY SYMMETRIC SOLUTIONS 



In the case of spherically symmetry the integration over angles in eq. (|2 . 16[) can be performed explicitly, yielding 
the one-dimensional non-linear integro-differential equation 



d 



2 7r rj 



dr' r' e* (r,) 



r — a\ 



where now 



Q = 47r / dr r e 



2 „Mr) 



(3.1) 



(3.2) 



Ea. l3.ll is to be supplemented with the boundary conditions of smooth joining (continuity of 4>{r) and <f>'(r) at r = l+a) 
with the external monopole solution <j>(r) = T}/r. 

Integration over the angles in eq. (|2.14p yields the non-linear integral equation 



(r) 



4 7T 77 



1 



1 



I 7 " - a \ 



max(r, a) J Q 



dr' r a e^" 1 + 0(1 - r - a) / dr' r' e^ r "> 

J r+a 



— / dr' r' 

2r J lr _ al 



r\2 



a . (r — r') 

— — r — r 

2 2a 



(3.3) 



6»(1 -r') 



which by itself determines 4>(r) for all values of r. In particular eq. (|3.3|) implies that any solution is finite in the 
origin, since 



Moreover, eq. (|3.ip or eq. (|3.3[) imply 



4tt 77 rl f a 



a Jo 



drr 2 e^ r) + / dr r e* (r) 



6'(0) = and <j)'(r) < for r > 



(3.4) 



(3.5) 



so that the density is monotonically decreasing away from the origin. The simplest way to prove that <fi'(0) < is to 
multiply both sides of eq. (|3.ip by r and integrate from to an arbitrary value. This yields 



<f>'(r) = 



2 7T 77 



dr r 



r +a 



dr" r" e^ r ") 



r' — a\ 
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which has a negative definite r.h.s. which vanishes when r — > as long as a > 0. 
The action for spherically symmetric solution becomes from eq. (|2.18|) . 

»M = ^ I r2 dr <K0 e 0(r) - log ^ . (3.6) 
y Jo 47r 

IV. SOLUTIONS REGULAR AS A -» 

The saddle point equation (|2.13[) admits solutions which are regular as a — > and reproduce in this limit the 
solutions of the standard Lane-Emden equation. We call them dilute solutions. These are known since longtime in 
the spherically symmetric case Q. To study such solutions we can set the cutoff A to zero from the start. Equivalently, 
one can take the limit a — > in eq. (|3.1[) with the assumption that <p(r) stays finite in the limit for any r, including 
the origin r = 0, and that </>'(()) = 0. 

Assuming spherical symmetry we get <fi(r) = rj/r for r > 1, while for r < 1 we have 

4>"{r) + - <j>'{r) + ^ e*« = , 0'(O) = , (4.1) 

with the boundary conditions 

0(1) = 77 , <t>'(l) = -r). (4.2) 

Clearly there exists one and only one solution for a given Q, so the classification of all solutions is equivalent to the 
determination of all the allowed values of Q. 
We set as usual @ 

4>{r) =0(O)+ X (Ar) , 0(0) = log ^ , X (0) = . (4.3) 
where xi z ) rnust satisfy, upon inserting eq. (|4.3[) in eq. (|4.1[) . 

^) + -x'W + e *»=0 , x'(0)=0 (4.4) 
and, from the boundary conditions eq. (|4.2[) at r = 1, 

log^+ X (A)=7 ? , A X '(A) = -?;. (4.5) 

These two relations fix the two parameters Q and A as functions of rj. In particular, we can rewrite <p(r) as 

(j>{r-ri) = V - X(KV)) + x(Hv)r) ■ (4.6) 
The density profile of eq. (|2.15[) reads in terms of x(A?"), 

p{r) = -f- e^ r ~> , (4.7) 

4 7T 1] 

so that the normalization factor can also be written Q = e v /p(l). 

Due to scale-invariant behaviour of the gravitational interaction, eq. ()4.4j) enjoys the following scale covariance 
property: if \{ z ) is a solution, then also Xa(z) = x( zea ) + 2a is a solution. Using x«( z ) rather than x( z ) is 
compensated by the shift A — » Ae _a . Notice that 0(r) is indeed defined as a scale transformation of x(z), with 
the scale parameter A (not uniquely) fixed by Ax' (A) = —rj. As a consequence of this scale invariance all physical 
quantities must be invariant under the simultaneous replacements A — > A er a and x(A) ~> x(A) + 2 a. 

Let us observe that, if the solutions of eq. (|4.1j) or (|4.4j) did not fulfill <fi'(0) — or x'(0) = 0, a delta function at 
the origin would appear in the right hand side and x( z ) — ~ A/ with A > 0, when z — > 0. This is the only possible 
singular behaviour and corresponds to a point particle with negative mass —A in the origin. We shall not consider 
such unphysical solutions, which are ruled out when eq. (|4.1[) is considered as the a — > limit of eq. (|3.ip . and stick to 

x'(o) - o. 
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Moreover, scale invariance allows to set also x(0) = 0. This choice completely fixes the scale and we have no more 
scale invariance left. One has from the second eq. (|4.5[) that x'i z ) < f° r au positive z, so that the density profile is 
monotonically decreasing with the distance. 

Since x( z ) 1S monotonically decreasing, the allowed values of Q [that is of p(l), see eqs. (12.151) and (j4.2|) ] are in 
one-to-one correspondence with to the roots of the relation A = X(rj) which inverts the second of eqs. (|4.5[) . Thus 
there is only one dilute solution for each value of A. However, A = X(rj) is a multiple valued function for a certain range 
of 77. We plot i] vs. log A in Fig. [T] One sees that there is a unique A for a given 77 only for A < 3.6358865 . . ., that is 
7] < 772 = 1.84273139 . . .. For 772 < V < Vc = 2.517551 . . . (r]c is the absolute maximum, located at Ac = 8.99311 . . ., 
of — Ax' (A) over < A < 00) the relation A = A (77) is indeed multivalued and as 77 approches the value 2 there are 
increasingly more solutions which accumulate near the purely logarithmic solution 




FIG. 1: rj as a function of log A according to eq. (|4.5[) for the dilute solutions. Notice the maximum of 77 at r\c = 2.51755 . . ., 
where cy diverges. The region beyond the maximum location, at log Ac = 2.196459..., is locally unstable as discussed in 
ref.0. cp and kt diverge before the maximum, at the point t\t = 2.43450 — 



0ooO) = lim 77 - x(A) + x(Ar) = 2-2 log r. 

A — >oo L 

which correspond to the singular density profile Pca(r) = (Awr 2 ) . In other words, when 77 is exactly 2 there exist 
an infinite set of solutions corresponding to the increasing sequence of values of A which satisfy Ax'(A) = —2. In Fig. 
[5] we plot in log-log scale the density profile p(r) [see eq. (|4.7[) ] vs. r for few terms of this sequence. 

We recall that, for any given 77 < 77c, only the unique dilute solution with A < Ac is locally stable. That is, 
any small fluctuation increases the action. At A = Ac a zero-mode appears in the linear fluctuation spectrum and 
cy, the specific heat at constant volume, diverges. For larger values of A negative modes show up, signalling local 
instability. Actually the locally stable dilute solutions describe a (locally) stable dilute phase, in the thcrmodnamic 
and mechanical terms, only for 77 < 777- = 2.43450... At A = A^ = 6.45071 . . . the isothermal compressibility diverges 
and the speed of sound inside the sphere becomes imaginary. The gas phase collapses at this point in the canonical 
ensemble as a consequence of the Jeans instability. This is confirmed by Monte Carlo simulations 5]. Hence, in the 
case 77 = 2 depicted in Fig. ^ only the solution with the smallest density contrast, corresponding to the smallest 
value A = 4.071496 . . ., is stable. The unstable solutions that exist for 772 < 77 < 777- should imply that the dilute phase 
is metastable for this range of 77, However, their lifetime is so huge [see ea. (l5.20p ] that these dilute solutions are in 
practice stable. 

In the stable dilute phase the particles are moderately clustered around the origin with a density that monotonically 
decreases with r. One sees that the density contrast between the center and the boundary 



pii) 



grows with A since xM decreases with A [see eq. 63}] and Fig. [U [3], Si]. 
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FIG. 2: Log-log plot of the density profile for some dilute solutions at r\ — 2. Only the smallest one at A = 4.071496 ... is 
locally stable. 



The action per particle of a dilute solution can be written entirely in terms of the density at the border Q 

"4-7T 



s(t?) = 3 - i] + log 



- 4ttp(1) = 3 - r, + log + X (A) - — e*W 



(4.8) 



and is multivalued for ?/ 2 ■ • • < V < Vc • ■ ■■ The solution with the smallest density contrast has the smallest action, in 
accordance with the linear stability analysis. We plot s(rj) as a function of log A and 77 in Fig. [31 
All physical quantities can be expressed for large TV in terms of the function 



„, . 7? ds 47T . , 



(4.9) 



We have for the energy, the pressure, the isothermal compressibility kt, the specific heats (cy and cp) and the speed 
of sound at the boundary Q 



NT Jy " V N 



S — So 

N 

Cy = 3 

vj 
T 



-3[l-/(ry)]- S (ry) + 



1 



3 [m - n f'(v) - 



o 



Cp = c v 



o 



m + 1 v f(v) 

im - V fir))? 



1 

N 

o 



(4.10) 



- + fa) + 7i V f'(v) + 0[- 



Here So stands for the entropy of the ideal gas. 

Notice that kt and cp, the specific heat at constant pressure, both diverge before 77c, at 77 = rjx = 2.43450 . . .0]. 
Moreover, the speed of sound squared, which stays regular at the boundary when 77 — > 77T, has a simple pole there 
when evaluated inside the isothermal sphere^]. Thus, when 77 > 77^ the speed of sound is purely imaginary inside the 
sphere implying that small fluctuations grow exponentially with time. 

We have here discussed the dilute solutions in the zero cutoff limit (A = 0). Including the short distance cutoff 
only adds small O(A) corrections to the dilute solutions. 
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FIG. 3: The action per particle s vs. log A (above) and vs. 77 (below) for the dilute solution. s(j]) is multivalued for rj > 772, 
with a peculiar self-similar structure, as evident from the inset. 



V. SOLUTIONS SINGULAR AS A -» 0. 



We present in this section new solutions for which the mean field <j>(r) develops singularities in the limit A — > 0. 
In case of spherical symmetry, when eq. (|3.ip is the saddle point equation, the only possible singularity is localized at 
r = 0. It is convenient then to 'blow up' the region near r = by setting [compare with eq. (|4.3[) ] 



-) = 0(0) + /i 2 £ (-) , 0(0) = log 9*. , £(0) = , so that p(r) 
\ a J inn 



A 2 



Attt] 4tt 77 

where we introduce the variable: 

/j = A a . 

We thus obtain for £(x), x = r / a, the following rewriting of the integral equation (|3.3[) 



^ 2 «(-) 



1 



max(x, 1) 
a; + 1 



I 2 (\x-l\)-h(m(x)) 



(--I) 2 



4.r 



h{m(x)) - h(\x - 1\) 



2x 



h{m{x))-I 2 {\x-l\) 



1 

4x 



h{m[x))-h{\x-\\) -/ 2 (l)+/i(l) 



where, 



m(x) = min f~,a; + 1^ , J„(a;) = ^ y™ exp[/j 2 



(5.1) 



(5.2) 



(5.3) 
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and we have used eq. (|3.4[) to write 0(0) in terms of £(x) as 



0(0) - /i 2 



h(l) 



(5.4) 



Notice that all dependence on A appears now through the variable [i = A a. 

By construction we have £'(0) = [see eq. (|3.5|) ]. £(0) = and £'(x) < 0, £(x) < for all x > 0. Notice that 77 does 
not enter the integral equation above, just as it did not enter the differential equation (|4.4|) for \ m the setup without 
short-distance regulator. Owing to the normalization condition eq. (|3.2p . 77 is computed once a solution is known, 
as a function of /j (and a), through 



77 = /J ah 



= fi a 



dxx 2 exp[/x 2 £ (x)] 



^From eq. (|3.1[) we may also derive the integro-differential form of the equation satisfied by £(2;), that is 



d 2 



/i(x + l)-Ii(|a:-l|) 



(5.5) 



(5.6) 



Eq. (|5.2j) or (|5.6[) can be solved numerically to high accuracy for any given /.i. We plot the corresponding 77 as a 
function of log A in Fig. [?J for better comparison with its behaviour in the dilute case, Fig. [TJ 




FIG. 4: ?7 as a function of log A according to eq. (|5.5[) . The black circle mark the values of A corresponding to the interpolating 
solutions depicted in fig. [S] Here, /i 2 = A 2 a 2 = C(log i). 



We find four different regimes according to the value of \i for a<Cl. 

• /i = 0(a). That is, A = 0(1). In this regime we reproduce the curve r\ = 77(A) (Fig. [T]) characteristic of the 
diluted phase of sec. IV plus small corrections of order a. 

• l-i = 0(1). Here A = O(^) and the damped oscillation pattern around 77 = 2, characteristic of the dilute 
solutions, gets disrupted: the oscillations regain larger and larger amplitude until /i 2 = A 2 a 2 ~ log -. 

• /1 2 = O(logi). Here, a sudden drop to very small 77 = rj m i n (a) ~ a takes place [see fig. 0]; then 77 rises again 
as [i a, for [i a large enough, and keep growing indefinitely (see fig. 2] and below). This regime corresponds to 
saddles interpolating between the dilute phase and a collapsed phase (see below). 
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• /i 2 = O(-). This regime describes collapsed configurations where all particles are concentrated in a region of 
size ~ a and provide the absolute minimum of the free energy. 

Notice that new maxima appear in addition to those present for a = (compare with Fig. [TJ). Those new maxima 
turn to be degenerate with the old ones up to small corrections of order a. In addition, 77 exhibits oscillations for 
/1 = O(l) which are similar to those present at a = and 77 — ► 2 but reversed and faster. We shall give below a 
qualitative explanation for this peculiar behavior. 

The new dependence of 77 on A described above implies that the number of solutions (that is the number of distinct 
A for any given 77) is finite for any 77 and that, most importantly, there appear new solutions. In particular, there 
are two new solutions for ?7 m in(a) < ?7 < 772 in addition to the unique solution regular as a — > and there is one new 
solution for 77 > 77c, the region not accessible to dilute solutions. These new solutions are singular in the limit a — » 0, 
as evident from Fig. [H where we plot 4>{r) vs. r for few values of a at a fixed value of 77 and A placed in the sudden 
drop of rj vs. log A (the black circles in Fig. IH scaling as /j 2 = A 2 a 2 ~ log ^). 



35 
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logio r 

FIG. 5: Profiles of the mean field 4>{r) for fi 2 = O(logi) (interpolating solutions), with comparison between the regularized 
and unregularized cases. The dashed black curve is the dilute solution (that is for a — 0) at the same value of 77, namely at the 
location of the magenta circle of Fig. [4] The other curves correspond to the black circles in fig. [4] 

In Fig. [5]we also compare the solutions of the integral equation (|3.3[) or (|5.2p to those of the differential equation 
(I4.1|) , with boundary conditions (|4.2[) , in which Q takes the value obtained by the integral method, upon using eqs. (|5.1|) 
and (|5.4p . We see that the agreement is very good down to r of order a, when the short-distance regularization becomes 
effective, while the solutions of the differential equation (]4.ip eventually blow to —00 as r — > 0. This agreement appears 
very natural from the integro-differential formulation of eq. (|3.1[) . where the short-distance regulator a can play a 
significant role only for r < a. 

Hence, the qualitative behaviour of the singular solutions can be understood already from the differential equation 
(|4.4p . provided one drops the initial conditions x(0) = x'(0) = characteristic of the dilute solutions. That is, using 
the relation x(//x) = /J 2 £(x) that follows from eqs. (|4.3[) and (|5.ip . we can solve the integral equation (|5.2p just in 
the interval < x < x\, with X\ > 2, to compute 

x(2/ J ) = M 2 e(x 1 ) and X '(2/x) = n^'( Xl ) , (5.7) 

at the boundary x = x%. These values for x(2 /i) and x'(2 /i), both negative by construction, can then be used as initial 
conditions at z = 2 fi to solve the differential equation (|4.4p in place of the conditions x(0) = x'(0) = 0, characteristic 
of the dilute solutions. The choice x\ > 2 is motivated by a direct numerical analysis but can be explained by the 
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FIG. 6: Profiles of the logarithm the density logp(r) for /j 2 = O(^) (collapsed solution). For the value of r\ given, we have 
A = 6012149.6972 .... This collpased profile should be compared to the dilute ones (stable and unstable) depicted in Fig. [2] 



a = 10" 



fact that fi diverges when a — > in the solutions which are singular as a — > 0. This implies, on one side, that the 
r.h.s. of eq. (|5.6[) is exponentially small in /j, 2 when xo > 2 [see also eq. (|5.3[) ] and, on the other side, that the initial 
conditions at z = 2 fi affect the function x{ z ) very strongly, causing the rapid drop for r > 2 a observed in Fig. [5l as 
well as the growth followed by the fall to — oo as r — > 0. It must be observed, in fact, that £(x) has a finite limit as 
(j, — > oo [see below eq. (|5.9p ]. so that to leading order and fj. ^> 1, using eq. (|5.7p . 

X (2 M )=C>( M 2 ) and X '(2 //) = 0( M ) . 

It follows that the function xi z ) fulfills approximately the free equation (zx)" = in the rapid drop for z > 2 fi as 
well as in the following plateau (to the right in Fig. [5]), that is 



X(z) 



Co 

z 



c = O( M 3 ), Cl = cV) 



(5.8) 



Where we used that z = A r = (X x. This behaviour is valid as long as x"( z ) e x<yZ \ that is z < C(/i exp[/x 2 /3]) 
since x = Oiv?) an d Co = C(/i 3 ). For larger values of z the function x(z) tends to the unique large distance fixed 
point given by the purely logarithmic solution log(2/z 2 ) and does this through the characteristic damped harmonic 
oscillations in logA0 responsible for the waving behaviour of rj in the dilute phase [see fig. [T]. Since the crossover 

to the logarithmic behaviour takes place for larger and larger values of z the smaller is a, because fi = O ^^log 

grows with a as well as the crossover point z — 0(/j, exp[/i 2 /3]). When x{ z ) is almost constant, the characteristic 
oscillations tend to have the same amplitude they have in the regular solutions. Moreover, this almost constant value 
of x{ z )i which is of order /i 2 , decreases faster in /j than the end point A = /i/a at which x( z ) is to be evaluated to 
give rj as in the second of eqs. (|4.5p . This explain why the oscillations of r\ are reversed and faster right before the 
sudden drop in Fig. [4j In other words, the oscillations in the left part of fig. [4] describe the approach to the limiting 
solution log(2/2; 2 ) [77 = 2] while one is getting off this limiting solution in the right part of fig. |3J This explains why 
the oscillations on the right get reversed compared with the oscillations on the left of fig. 0] 

In the next two subsections we will provide more details on all these matters and further analytical insights. 



A. Collapsed phase 

This is the case when \i ~ rj/a as a — > and corresponds to the growing branch on the right of Fig. [U Notice 
that we have collapsed solutions for arbitrarily large values of 77 contrary to the dilute solutions which only exist for 
f] < r) C - 

Since is negative and monotonically decreasing for x > 0, the quantity exp[/j£(a;)] is exponentially small except 
in the interval < x < xq where £(x) would vanish in the limit. We use now this property to get an approximate 
but analytic singular solution. That is, the dominant contribution in the r.h.s. of eq. (|5.2p is obtained when all 
integrations are restricted to the interval < x < xq. One easily realizes that, by consistency, xq = 1/2. Then one 
finds 



I n {x) 



1 



n+l 



x < 1/2 
x > 1/2 
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and in particular using eas. (|5.4p and (|5.5 



0(0) = /i 2 



a 



h(l) 



Finally, the leading form for 0(r) reads 



24?? rr 

0(o +— d- 

a \a 



EL 

24 



1 
a 

V V_ (t_ _ 1 

a 2r\a 2 

r 



/j ah - 
1 a 



/^ 2 a 
~2T 



0(0) = 2 



r < a/2 

a/2 < r < 3a/2 
r > 3a/2 . 



(5.9) 



Notice that 0(r) as well as its first and second derivatives are continuous at r = a/2 and r = 3 a/2. For r > 3a/ 2, the 
effective gravitational field has the free form proper of a locally vanishing mass density. 

Higher corrections can be obtained by iteration over eq. (|5.2[) . That is, using the leading form of £ (x) that can be 
read out of in eq. (|5.9p to evaluate to next-to-leading order, as /1 — > 00, the integrals I n (x) in the r.h.s. of eq. (|5.2p . 
This procedure can then be repeated to obtain further corrections. 

Eq. (|5.9[) provides the spherically symmetric, leading order singular solution of the mean field equation (|2.13[) where 
the particles are densely concentrated around one point in a region of size of the order of the cutoff A. In fact, to 
leading order in a the density reads 



(i 



while to next-to-leading order we have 



p(r) 



p(r) 



< a/2 



7r a J 

. , r > a/2 , 



■ exp 



(5.10) 



where 0(r) has the leading order form of eq. (|5.9[) . This corresponds to the situation when the gravitational attraction 
completely overcomes the kinetic energy and all particles fall towards the origin. Of course, such singular solution 
only exist mathematically for non-zero values of the cutoff A. 

Let us now compute the action for the solution eq. (|5.9p from eq. (|2.1ip . The regularized laplacian can be computed 
from the equation of motion (|2.13p with the leading order result for a — * 0, 

d 3 r 0(r) V 2 0(r) = -^LLH 4^ f r ? dr 0(r) e^ r) = -4tt 77 0(0) = . 

Q Jo a 



\r\<l 

The second term is just 



logQ = 0(O)-log 



A 



Att 77 a 



log : 



so that all together, to leading order we obtain for the free energy 

Gm 2 (3N Gm 2 P N NTij 8 
F - F Q = — — + NT s(r), a) = — — — + NT log — 



2.4 



2 A 2 a 

Gm 2 N(N - 1) 
2A 



NT log 



(5.11) 



(A/2)* 



The first term is just the potential energy of N particles clustered in a small sphere of radius A/2, where the regularized 
gravitational interaction is the same for all ^N(N — 1) particle pairs. The second term is T times the entropy loss 
in the collapse. As one could expect, this free energy is large and negative, unbounded from below as A — * 0. In 
particular, this free energy is well below the free energy of the dilute solution eq. (|4.8p and Fig. [3] for all values 77 > 0. 
That is, the singular solution eq. (|5.9p provides the absolute minimum for the free energy eq. (|2.1ip for any 77 > 0. 
The dilute solution eas. (|4.4p and (|4.2p is only a local minimum. In the next subsection we shall address the issue of 
the metastability of the dilute phase. 
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Let us now compute some physical quantities characterizing the collapsed phase. 
The equation of state has the mean field form 



PV 
~NT 



47T 



(5.12) 



in both the dilute and collapsed phase. In the former case the r.h.s. of eq. (|5.12[) coincides with the function f(rj). In 
the latter we obtain instead, from eqs. (|5.10p and (|5.9|) . to next -to-leading order 



PV 

iVT 



exp 



V 
a 



(5.13) 



Thus, the pressure at the boundary is exponentially small in the collapsed phase. The complete plot of PV/NT vs. 
r\ is given in Fig. [7] One can see how the well known [3J] spiraling towards (rj — 2, f (rj) — 1/3), characteristic of the 
dilute solutions, is now limited to the case when jj, is of order a . When [i reaches values of order 1 the inward spiral 
comes to a halt and then winds back until < C(log -). For larger values of /x, that is /i ~ T]/a as in the collapsed 
phase, PV/NT drops to values exponentially small in ~. Analogous plots are obtained for other types of cutoff in 
refs. [|. 




l x = 0{l) 



2.5 



FIG. 7: Complete plot of PV/NT vs. r\ when a = 10~ 4 (blue curve). When \i is of order a the blue curve cannot be distinguished 
from the corresponding curve of the case a = 0, here depicted in red (dilute solutions). When p is of order 1 the blue curve 
(a > 0) and the red curve separate. The blue curve then traces the red spiral backwards at a distance of order a (log a) 2 as long 
as /i 2 < C(log -) (interpolating solutions). Finally the blue curve drops to zero when /i 2 is of order i (collapsed solutions). 

For other physically relevant quantities we similarly find 



CV= 2 



c P = - + 0(a 2 ) 



vj 1 8 

Tf = — = -T exp 

1 kt a 



S — *5o 
N 



3 loga + O(a ) . 



(5.14) 



B. Interpolating solution: metastability 

The dilute solutions studied in sec. IV are metastable because the collapsed solution provides the absolute minimum 
for free energy [eq. (|5.1ip ] . The collapse transition from the dilute phase to the collapsed phase can take place through 
an interpolating solution in an analogous way to the bubbles in liquid-gas first order phase transitions [101 ] . 

Such interpolating solutions are the saddle points of eq. (|2.13|) , which correspond to the sudden drop of ry as a function 
of log A (see Fig. @| located approximately at A 2 et 2 ~ log -. In this narrow drop region we have 0(a) < r) < r)c, so 
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the interpolating solutions are allowed only in this interval of 77. Examples of the interpolating solutions are given in 

% El 

According to the qualitative discussion below eq. (|5.8j) , the location of the 77 drop can be obtained approximately as 
the values of A at which the function x(A) nas the crossover from the free solution x(A) — co/A + ci to the logarithmic 
solution log(2/A 2 ). That is for A 2 a 2 ~ log i, since c = 0(A 3 a 3 ) and ci = 0(A 2 a 2 ), see eq. I|5.8J) . Then the scaling as 
a — > of p(0), the density in the origin, can be obtained from eq. (|5 . 1 [) and the numerical analysis as 



P(0) 



4tt rj 



o(a- 2 io g i; 



because 77 = 0(1) in this regime. 

On the other hand, from Fig. [5] one sees that the value of the density for r ^> a and in particular the value at the 
border p(l) differ by order a w.r.t. the dilute phase. The equation of state (|5.12|) . which holds for any value of A, 
then implies that the pressure in the dilute and interpolating solution, for the same value of 77, differ only by order a, 
as evident also from Fig. \7\ 

In practice, one can regard an interpolating configuration as a special spherical fluctuation of the dilute phase 
in which a fraction of particles of order a log ^ is almost uniformely removed from everywhere and concentrated in 
a region of size a around the origin. In fact, this result has an heuristic explanation by a simple energy-entropy 
argument: if we denote by Ap the variation of the density in a region of size a around the origin, we then have 



77 /AN\ 2 AJV, 1 AN , A 

As + 3 log - , = a 3 Ap 

2a \ N J N B a N H 



(5.15) 



since the first term is the gain in potential energy per particle and the second the loss of entropy per particle due to 
the concentration. We may neglect the effect due to the rest of the particles if we assume that the fraction of 
particle moved around, would vanish as a — > 0. Minimizing eq. (|5.15[) w.r.t. Ap yields 



Ap 



3 1 1 

— a lo S _ 

77 a z a 



and As = 



9a 
2^ 



log 



(5.16) 



which is indeed consistent with the requirement that a 3 Ap vanishes as a — > and rather accurately approximate our 
numerical results. 

Notice that the minimized As is to be identified with the difference, between interpolating and dilute solution, of 
the action per particle, so that 



As = 



^interpolating ^dilute 



2r) 



a (log a) 2 



(5.17) 



This result is in very good agreement with the numerical values obtained by solving eqs. (|3.3[) and 

According to the interpretation of the interpolating solutions as saddle points through which the dilute phase may 
decay to the collapsed one, the lifetime of the dilute phase is given by 



f jDetj 
V Dct 



1/2 



NAs 



(5.18) 



where Det stands for the determinant of small fluctuations around the interpolating solution and Deto for that around 
the dilute solution. 

This ratio of determinants can be estimated as follows: in the interpolating profile there should be only one negative 
mode necessarily localized in a region of size a around the origin, where the density is very large, since —p(r) plays 
the role of Schroedinger potential in the eigenvalue equation for the small fluctuations. Eq. (|5. then implies that 
-i? log i provides the scale of the potential felt by this negative mode. Hence the corresponding negative eigenvalue 
should be of the same order, while the rest of the spectrum should be almost unaffected, since the two density profiles 
differ only by order a for r > 2 a. Therefore 



I Det I 
Det 



1 , 1 

-3. lo ST 



(5.19) 



We obtain from eqs. (|5.17| - (|5.19p to leading order the following lifetime for the dilute phase of the self-gravitating 
gas: 



1 



1 



a V a A 



log^ exp 



9 AT 



2Gm 2 



(log 



(5.20) 
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One can see that the lifetime becomes infinitely long in the zero cutoff limit as well as when N — » oo at fixed cutoff 
[recall that R ~ N in the dilute limit of eq. (|1.1[) ]. 

We want to stress that the lifetime of the dilute solutions is extremely long for large N and small cutoff A. This is 
due to the fact that although the collapsed solution has a much lower energy than the dilute solution, in the latter 
there is a huge entropy barrier against the gathering, in a small domain of size ~ A, of a number particles large 
enough to start the collapse. A qualitatively similar conclusion is reached in ref. [TTj |. 

We shall come back for a more detailed study of this problem in a subsequent paper. 
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